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Abstract 

We report the results of a study of neutral B meson decays to the D^tt^tt" final state, where the 
is fully reconstructed. The results are obtained from an event sample containing 388 million 
BB-meson pairs collected in the Belle experiment at the KEKB e'^e~ collider. The total branching 
fraction of the three-body decay D^Tr+Tr") = (8.4±0.4(stat.) ±0.8(syst.)) x 10"^ has been 

measured. The intermediate resonant structure of these three-body decays has been studied. From 
a Dalitz plot analysis we have obtained the product of the branching fractions for and D^~^ 
production: B{B^ DI+-K-) x B{Dl+ D^tt+) = (2.15±0.17(stat.) ±0.29(syst.) ±0.12(mod.)) x 
10-4, and B{B^ D*+n-) x 3(0*+ D^Tr+) = (0.60 ± 0.13(stat.) ± 0.15(syst.) ± 0.22(mod.)) x 
10"4. This is the first observation of the S° DQ'^-7r~ decay. The B^ D'^p^ and -D°/2 branching 
fractions arc measured to be: B{B^ D^p^) = (3.19±0.20(stat.) ±0.24(syst.)±0.38(mod.)) X lO^^, 
and B{B^ D°/2) = (1.20 ± 0.18(stat.) ± 0.21(syst.) ± 0.32(mod.)) x lO"''. 

PACS numbers: 13.25.Hw, 14.40Lb, 14.40.Nd 
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INTRODUCTION 

The decay D^tt^tt^ includes intermediate states Z}**^7r^, where D**'s are P-wave 

excitations of states containing one charmed and one hght (g = u, d) quark that decay to the 
D°7r+ final state. Fi gure 1 shows the spectrum and the allowed transitions of cg-meson states. 
In the heavy-quark limit, the c-quark spin Sc decouples from the other degrees of freedom, 
and the total angular momentum of the light quark jg = L + Sg is a good quantum number. 
Four P-wave states with the quantum numbers 0^{jg = 1/2), l~^{jg = 1/2), l~^{jq = 3/2) 
and 2~^{jg = 3/2) are expected; these are usually labeled as Dq, D[, Di and D^, respectively. 
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FIG. 1: Spectrum of cg-meson excitations. The lines indicate possible one-pion transitions. The 
Dq, D'l mesons are broad which is indicated by shaded areas. 

The two jg = 3/2 states have narrow widths of about 20-40 MeV and are well estab- 
lished [1-11]. The measured masses agree with model predictions [12-15]. The remaining 
jg = 1/2 states are expected to be broad and decay via S'-waves. The B Dnir decay 
process provides a way to study D** production. Angular analysis of the decay products 
can be used to determine D** meson quantum numbers. These results also provide a test of 
Heavy Quark Effective Theory (HQET) and QCD sum rules [16, 17]. 

A study of neutral D**^ production in charged P-decays has been recently reported by 
Belle [18], where four D** states are observed and the production rates of the broad {j = 1/2) 
states are found to be of the same order-of-magnitude as those for the narrow (j = 3/2) 
states. This paper describes an analysis of the P° P^tt+tt^ decay that is performed in a 
manner similar to that of the previous Belle analysis of the B^ P**+7r^ decay [19]. The 
results presented here supersede those of Ref. [19]. 

The neutral P decay to P**7r is described by the tree diagram only as shown in Fig. 2(a) 
while for the charged P decay to D**tt, the amplitude receives contributions from both 
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tree and color-suppressed diagrams as shown in Fig. 2(b,c). D** tree-diagram production 




FIG. 2: Quark-line diagrams for neutral (a) and charged (b and c) B decays. 



amplitudes are described by the Isgur-Wise functions ri/2 and r3/2. According to the QCD 
sum rule [16, 17], ti/2 ^ and one would expect suppression of decays to the broad state. 
The observation that the production rates of the broad (j = l/2)-states are comparable 
with those of the narrow {j = 3/2)-states indicates either a large contribution of the color- 
suppressed diagram or the violation of the sum rule above. Measurement of the decay rates 
of the neutral B allows one to test the contribution of the tree-diagram only and also test 
the QCD sum rule. 

In this analysis, the final state contains two pions of opposite sign, and these can orig- 
inate from resonant states such as the p°, /o, /2, etc. While the possible presence of vrvr 
resonant structures complicates the analysis, it can also provide valuable information about 
the mechanism of these decays. 



THE BELLE DETECTOR 

The Belle detector [20] is a large-solid-angle magnetic spectrometer that consists of a 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) for charged parti- 
cle tracking and specific ionization measurement (dE/dx), an array of aerogel threshold 
Cerenkov counters (ACC), time-of- flight scintillation counters (TOF), and an array of 8736 
CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux return located outside the 
coil is instrumented to detect K^^ mesons and identify muons (KLM). We use a GEANT- 
based Monte Carlo (MC) simulation to model the response of the detector and determine 
its acceptance [21]. 

Separation of kaons and pions is accomplished by combining the responses of the ACC 
and the TOF with dE/dx measurements in the CDC to form a likelihood C{h) where h = (vr) 
or {K). Charged particles are identified as pions or kaons using the likelihood ratio (7^): 

A more detailed description of the Belle particle identification can be found in Ref. [22]. 
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EVENT SELECTION 

A data sample of 357 fb~^ (388 million BB events) collected at the T(45') resonance is 
used in this analysis. Candidate S° £)°7r+7r~ events are selected, where the mesons 
decay via the — > K~7t~^ mode. The signal-to-noise ratios for other decay modes 
are found to be significantly lower and, therefore, are not used. (The inclusion of charge 
conjugate states is implied by default throughout this paper.) 

Charged tracks are selected with requirements based on the average hit residuals and 
impact parameters relative to the interaction point. We require that the polar angle of each 
track be in the angular range of 17° — 150° and that the track transverse momentum be 
greater than 50 MeV/c for kaons and greater than 25 MeV/c for pions. 

Charged kaon candidates are identified by the requirement Tl^K) > 0.6, which has an 
efficiency of 90% and a pion misidentification probability of approximately 10%. For pion 
candidates we require 7?-(7r) > 0.2. Kaon and pion candidates are rejected if the track is 
positively identified as an electron. 

Candidate mesons are K't:^ combinations with an invariant mass within ±12 MeV/c^ 
of the nominal mass, which corresponds to ~2.5 ok-k- We reject D° candidates that, when 
combined with any 7r° in the event, has a value of Md^^ — Mjjo that is within ±2.5 MeV/c^ 
of the nominal D*^-D^ mass difference. 

B meson candidates arc identified by their center-of-mass (cm.) energy difference Ai? = 

(JZiEi) — E]^, and the beam-constrained mass Mbc = \J E^ — {Y,iPi)'^i where E\^ = v^/2 is 
the beam energy in the T(45') cm. frame, and Pi and E^ are the cm. three-momenta and 
energies, respectively of the B meson candidate decay products. We select events satisfying 
Mbc > 5.25 GeV/c^ and \AE\ < 0.10 GcV. 

To suppress the large continuum background (e+e^ qq, where q = u,d, s, c), topological 
variables are used. Since the produced B mesons are almost at rest in the cm. frame, the 
signal event shapes tend to be isotropic while continuum qq events tend to have a two- 
jet structure. We use the angle between the thrust axis of the B candidate and that of 
the rest of the event (Othrust) to discriminate between these two cases. The distribution 
of I cos Othrust I is strongly peaked near | cosGthmstl = 1 for events and is nearly fiat for 
T(45') BB events. We require |cosGthrust| < 0.8, which eliminates about 83% of the 
continuum background while retaining about 80% of signal events. 

There arc events for which two or more track combinations pass all the selection criteria. 
According to MC simulation, this occurs primarily because of the misreconstruction of the 
low momentum pion from D** Dtt decays. To avoid multiple entries, the combination 
that has the minimum difference of z coordinates at the interaction point, l^^ri ~ -2712!) 
the tracks corresponding to the pions from B — >■ Di^i^n^ are selected [24]. This selection 
also suppresses combinations that include pions from decays. In the case of multiple 
D —>■ Kt: combinations, the one with the invariant mass closest to the mass is selected. 

^ D^TT+TT- BRANCHING FRACTION 

The D^7r'^Tr~ final state, together with three-body and quasi-two-body contributions, 
includes the two-body B^ — > D*~^7r~ decay followed by the decay D*~^ D^t:^ . We obtain 
the branching fraction of the three-body decay B^ D^vr+vr" excluding the contribution 
of B^ — > D*^Ti:~ . Using the Mb^k — Md mass difference, we subdivide the total sample in 
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to two subsamplcs as follows. Events that have a Dvr combination with Mdt^ — Mr, within 
3 MeV/c' (~ 6a) of the nominal D*+ - mass difference are denoted below as sample (2); 
the rest of the Dtttt events are denoted as sample (1). Sample (2) is used to crosscheck our 
procedures. 

The Afbc and AE" distributions for S° 7r"'"7r events arc shown in Fig. 3. The 

distributions are plotted for events that satisfy the selection criteria for the other variable: 
\AE\ < 25 MeV and |Mbc — ?«b| < 7 MeV/c^ for the and AE histograms, respectively, 
where me is the nominal B mass. Distinct signals are evident in all of the distributions. 

The background shape is obtained from generic MC data samples that include B^B~ 
(BC) and B^B^ (BN), continuum charm production (CC) and continuum with light quarks 
(UDS), each corresponding to approximately twice the luminosity of the experimental data. 
The Dir and tttt invariant mass distributions are different for different MC samples. The 
branching fractions used in the generic MC are measured with some experimental uncertainty 
and may not reproduce the experimental data. To improve the quality of the MC spectra, 
relative weights of these four components are determined from a fit to a two-dimensional 
Dalitz plot distribution for events in the AE sideband shown in Fig. 4. The fitting function 
represents the sum of the four two-dimensional histograms with floating weights. Each 
histogram is determined from its respective MC sample. The weights obtained for the four 
components are: a^c = 1-10 ± 0.07, gbn = 1-37 ± 0.22, acc = 0.52 ± 0.12, auns = 
0.92 ± 0.22. The AE background shape is described as Fbg{AE) = J2iO'iF'i{AE), where 
Fi{AE) is the AE distribution of the i-th component obtained from the MC sample. 

The signal yield is obtained by fitting the AE distribution to the sum of two Gaussians 
with the same mean value to describe the signal, plus the above-described background 
function Fhg{AE). The width of the broader Gaussian and the relative normalization of 
the two Gaussians are fixed to the values obtained from a MC simulation; the signal and 
background normalization as well as the width of the narrow Gaussian are left as free 
parameters. 

The fitted signal yields arc 2909 ±115 events and 4202 ± 67 events for samples (1) 
and (2), respectively. The reconstruction efficiencies (23.4 ± 0.4)% and (19.0 ± 0.4)% are 
determined from a MC simulation that uses a Dalitz plot distribution that is generated 
according to the model described in the next section. Taking into account B{D^ — > K~7r~^) — 
(3.80 ± 0.07)% [23], we obtain the following branching fraction: 

B{B^ D\+7r-) = (8.4 ± 0.4 ± 0.8) x 10-^ 

where the first error is statistical and second error is systematic. Various contributions 
to the systematic error are listed in Table I for both samples. They include tracking effi- 
ciency, particle identification efficiency, limited MC statistics, and background uncertainty. 
The background uncertainty is obtained_by varying the relative weights within their errors. 
The contribution of the non-resonant B^ — > K'tt'^tt'tt'^ is estimated using the Ktt mass 
sidebands of the D mass region and is negligible. 

The value of B{B^ — > D^tt+tt") improves and supersedes previous Belle result i3(i?° —>■ 
D^TT+TT-) = (8.0 ± 1.6) X 10-^ [25]. The value of the branching fraction B(S° ^ D*+7r-) is 
obtained usmg sample (2) and the PDG value B{D*+ D^n^) = (67.7 ± 0.5)% [23]. The 
resuh is i3(i?° ^_D*+n-) = (2.22 ± 0.04 ± 0.19) x 10"^ which is somewhere lower than the 
CLEG result B{B^ D*+7t-) = (2.81 ± 0.25) x lO'^ [26]. 
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FIG. 3: Afbc and AS distributions for 5° vr+vr events. Sample (1) distributions are shown 

in (a) and (c); those for sample (2) are shown in (b) and (d). 

B — > Dtttt Dalitz plot analysis 

For a three-body decay of a spin zero particle, two variables are required to describe the 
decay kinematics; we use the D^tt^ and tt+tt" invariant masses squared, and M^^, 

respectively. 

To analyze the dynamics of -B ^ Dnn decays, sample (1) events with AE and Mbc within 
the signal region {{AE + ^(Mbc — rnB))/<JAE)'^ + ((Mjc — ^^b) I ctm^^Y < 4 are selected. The 
parameters ctab = 11 MeV, (Tm^^ = 2.7 MeV/c^ and k = 0.9 are determined from a fit to 
experimental data; the coefficient n accounts for the correlation between Mbc and AE. 
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TABLE I: Tlie systematic uncertainties for tfie B D tt^tt branching fraction measurement. 

To test and correct the shape of the background, we use events from the AE sidebands, 
which are defined as: {{AE ± 65MeV + ^(Mbc - mB))/(yAEf + ((Mbc - mB)/ciM^,f < 4. 
Figure 4 shows the signal and sideband regions in the M\^c-AE plane. 

o 

H0.05 




-0.05 



-0.1 

5.25 5.26 5.27 5.28 5.29 

M^^ (GeV/c') 

FIG. 4: The experimental distribution in the (Mbc-AE^) plane. The ellipses show the signal (1) 
and sideband regions (2 & 3). 

The Dt[ and vrvr mass distributions for the signal and sideband events (sample 1) are shown 
in Fig. 5. In the Dtt mass distribution a narrow peak corresponding to D2 is evident. The 
TTTT mass distribution shows a peak corresponding to the p meson as well as a structure at the 
1.2 — 1.3 GeV/c^ mass region that is presumably due to /o(1370) and /2(1270) production. 

The M})^ and M^^ Dalitz plot distributions for the signal and sideband regions are shown 
in Fig. 6. The Dalitz plot boundary is fixed by the decay kinematics and the masses of the 
daughter particles. In order to have the same Dalitz plot boundary for both signal and 
sideband event samples, fits where the Kn mass is constrained to M/5 and Dtttt mass to 
are performed. The mass- constrained fits also slightly improve the accuracy of M|,^ and 
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FIG. 5: Dtt (a) and tttt (b) mass distributions for sample (1) events with helicity angles cos 9h > 0. 
The points with error bars correspond to the signal region events, the hatched region indicates 
the background obtained from generic MC events normalized to the sideband data. The open 
histogram shows the fitted function after efficiency correction. 



To extract the amplitudes and phases of different intermediate states, an unbinned fit 
to the Dalitz plot is performed using the method described in Ref. [18]. The event density 
function in the Dalitz plot includes both the signal and background functions. 

The backgrounds in the Dalitz plot are mostly combinatorial and have neither resonant 
structure (Fig. 6b) nor specific helicity behavior. The background shape is obtained from 
an unbinned fit to AE sideband events using the weights described above. The background 
Dalitz plot density is modeled by a smooth two-dimensional function. The number of back- 
ground events in the signal region is scaled according to the relative areas of the signal and 
sideband regions. 

There is no general method to describe a three-body amplitude. In this paper we rep- 
resent the Dtttt amplitude as the sum of Breit-Wigner function contributions for different 
intermediate two-body states. Such an approach is not exact because it is neither analytic 
nor unitary and does not take into account a complete description of the final state inter- 
actions. Nevertheless, the sum of Breit-Wigner functions describes the main features of the 
amplitude behavior and allows one to find and distinguish the contributions of two-body 
intermediate states, their interference, and the effective parameters of these states. We 
followed the same approach in the analysis of charged B decays [18]. 

In the D'^vr+vr" final state, a combination of the D°-meson and a pion can form a vector 
meson D*~^, a tensor meson Dg"*" ^ scalar state -Dg"*"; the axial vector mesons Df and D'^ 
cannot decay to two pseudoscalars because of angular momentum and parity conservation. 
The region of D^tt~^ invariant mass that corresponds to the D*~^ is excluded from the fit by 
requirmg \MDn - Md' I > 3 MeV/c^. However, in B -meson decay, a virtual D*~^ (referred to 
as D*) can be produced off-shell with Md-,^ above the D^7r~^ production threshold and such 
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FIG. 6: The Dalitz plot for (a) signal events; (b) sideband events. 



a process can contribute to the amplitude. Another virtual hadron that can be produced 
in this combination is B*~ (referred to as B*): B BI-k and B* —>■ Dn. For the mass of 
B*^ as well as the mass and width of the D*~^, we use the PDG values [23]; the widths of 
B*~ are calculated from the width of the D*~ in the HQET approach. To describe the vrvr 
system we include p, u, /2(1270), and three scalar mesons /o(600), /o(1370) and /o(980). 
The masses and widths of the p, u and /2(1270) mesons are fixed at their PDG values; the 
parameters of the scalar mesons are taken from the published papers on the /o(600) [28], 
/o(1370) [29] and /o(980) [29] . 

The contributions from the intermediate states listed above are included in the signal- 
event density qf)) parameterization as a coherent sum of the corresponding amplitudes 
together with a possible constant amplitude (ops). The phases of the amplitudes are defined 
relative to 

Siq',ql) = \aj,.A^*^{q\ql)+ao'/HA^*^{q\ql)+an,e^K 

+ape'^^{AP{q^ q^) + r^^.e'^- A-{q\ qf)) + aj.e'^^f^^*^^ Af^{q\ qf) 

+af,^u70)e'^'^'<^^''''^-''^^^Af'^^'^^^^ (1) 

where q^ = M|,^ and qf = M^^. The relative amplitude and phase of the u meson are 
expressed via those of the p meson. The relative phase is taken from p-u interference 
measurements [30], and the relative amplitude (r^^-p = aui/dp) is recalculated using that 
value. Assuming that the p and uo mesons produced in B^ decay emerge from the dd pair, 
the relative amplitude is expected to satisfy a relation rt^^p(^B) = —'^f^-pi'^)- 

We use the approach described in [18], where each resonance is described by a relativistic 
Breit-Wigner function with a q^ dependent width and an angular dependence that corre- 
sponds to the spin and parity of the intermediate- and final-state particles. The p meson 
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amplitude is described by the Gounaris-Sakurai parameterization [31]. We take into account 
transition form factors for hadron transitions using the Blatt-Weisskopf parameterization [32] 
with a hadron scale r=1.6 (GeV/c)~^. 

The variation of the detection efficiency over the Dahtz plot is taken into account by 
the minimization procedure. The efficiency dependence enters the likelihood function only 
through the normalization term. The normalization is obtained based on a large MC sam- 
ple generated uniformly over the Dalitz plane, processed with the same selection criteria as 
the data and multiplied with the model used to fit the data. The detector resolution for 
the invariant mass of the Dn{7nr) combination is about 2.5 (3.5) MeV/c^, which is much 
smaller than the narrowest peak width of 30-40 MeV/c^. Hence convolution of the de- 
scribed parametrization with the resolution is not necessary. The mass and width of the 
broad {Dn) resonance, M^*+ = 2308MeV/c^ r° + = 276MeV/c2 are taken from our D**^ 

^0 

measurement [18]. 

Table II gives the fit results for different models. The contributions of different states are 
characterized by their fractions, which are defined as: 



a 



I\Y-kake"t'^MQ)\'dQ' 

where Ai{Q) is the corresponding amplitude, and and (pi are the amplitude coefficients 
and phases obtained from the fit. The integration is performed over all available phase 
space characterized by the multidimensional vector Q (for decay to 3 spinless particles, 
dQ = dq'^dqf), and i is one of the intermediate states: D2, D^, p, /2, /o, D*, B* or 
the constant term a^g. The sum of the individual fractions fi exceeds unity for our case 
because of destructive interference. The product of the branching fractions of the B meson 
is expressed via the fraction ff 

^B^iTv^i-*DTv — _ 1 (3) 

where A^sig is the efficiency corrected number of the reconstructed Dtttt events and N^b is 
the number of BB pairs produced. 

Table II contains information on the likelihood change relative to the main set, and 
values obtained from four histograms: projections of M^^^ and M^^ for negative and positive 
helicities of the Dn and tttt systems, respectively[34]. 

The fit gives a statistically significant contribution from off-shell D*7r production; the 
addition of the off-shell S* amplitude does not improve the likelihood value significantly. 
The inclusion of the three-particle phase space term improves the likelihood value, however 
there is no reason to expect a constant amplitude with no momentum dependence over 
such a wide range of final particle momenta. Table III shows that the hkelihood changes 
significantly when the broad resonance is removed or treated as either a vector or a tensor. 
The change of likelihood — 21n£/£o = 51 for 2 additional degrees of freedom (amplitude 
and phase of Dq) corresponds to a significance of 6.8 a [33]. 

The branching fractions of D2 and Dq remain constant within errors for different models. 
The set of states used for the final results are D^, Dq, D*. p, /2 and the three above-listed 
/o's listed above, corresponding to column 1 in Table II. 

The values of the 1)2^ resonance mass and width obtained from the fit are: 



M^*+ = (2465.7 ± 1.8 ± 0.81^;?) MeV/c^ r^*+ = (49.7 ± 3.8 ± 4.1 ± 4.9) MeV, 
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Hp.(10-4) 











0.33 ±0.18 







629/603 


680/605 


632/601 


618/601 


659/605 


CL (%) 


23 


1.8 


18 


31 


6.3 



TABLE II: The fit results for different sets of amplitudes. Option 1 is used as the main set of 
amphtudes for the final results. 
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no broad state 


0+ 


1" 


2+ 


-2 In/: /A 


51 





28 


27 


CL, (%) 


659/605 
6.3 


629/603 
22 


652/603 
8.1 


640/603 
14 



TABLE III: Changes in likelihood for different quantum number assignments for the broad reso- 
nance. 

where the third error is model uncertainty. These parameters are consistent with but more 
precise than previous measurements performed by CLEO Mjj*q — (2463± 3± 3) MeV/c^ [8] 
and FOCUS To* = (34.1 ± 6.5 ± 4.2) MeV [27]. 

The product of the branching fraction for production obtained from the fit is: 

B{B^ D*+T^-) X B{D*+ D\+) = (2.15 ± 0.17 ± 0.29 ± 0.12) x 10'^ 

where the three errors are statistical, systematic, and a mo del- dependent error, respectively. 
We observe the production of the broad scalar Dq^ state with the product branching fraction, 

B{B° D*+7r-) X B{D*+ D\+) = (0.60 ± 0.13 ± 0.15 ± 0.22) x 10"^ 

This is the first observation of this decay (the interpretation of the neutral partner of this 
state is still a subject of theoretical discussion [35]). The relative phase of the amplitude 
is 

00 = 3.00 ±0.13 ±0.10 ±0.43. 

The Dtt helicity angle distributions for Mdtt regions corresponding to the D2 and are 
shown in Fig. 7(a) and (b), respectively, together with the efficiency-corrected fitting func- 
tion. The histogram in the region of the meson clearly indicates a D-wave. The distri- 
butions in the other regions show reasonable agreement between the fitting function and the 

data. 

The uncertainty of the background is one of the main sources of systematic error. This is 
estimated by comparing the fit results for the case when the background shape is taken sepa- 
rately from the lower or upper AE sidebands. The fit is also performed with more restrictive 
and looser cuts on AE, Afbc and AM^ that change the signal-to-noise ratio by factors of 
about two. The results obtained are consistent with each other and the maximum difference 
is taken as an additional estimate of the systematic uncertainty. The systematic errors on 
the Bi measurements (Eq. (3)) for the individual intermediate states include uncertainties 
in track reconstruction and PID efficiency, as well as the error in the K~ti^ absolute 

branching fraction. The model uncertainties are estimated by comparing fit results for the 
case of different models and for values of the parameter r of the transition form factor [32] 
fi-om to 3 (GeV/c)-^ 

The TTTT helicity angle distributions for M^^r ranges corresponding to the p, /2 and the 
region below the p, where the broad resonance dominates, arc shown in Fig. 8. For the 
positive helicity region, where the Dti contribution is suppressed, a clear P-wave structure 
for the p and D-wave structure for the /2 is observed. The scalar component parameters 
cannot be determined from the fit. This process can also have contributions from non- 
resonant background. The product branching fraction for the /2 is 

B{B^ D^f2)B{f2 TT+TT") = (0.68 ± 0.10 ± 0.12 ± 0.18) x 10"^ 
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FIG. 7: Dtt helicity angle distributions for data (points) and MC (histogram). The hatched 
distribution shows the background from the AE sideband region with the appropriate normal- 
ization, (a) corresponds to the D2 region {Md-j^ — 2.46| < O.lGeV/c^; (b) the Dq region 
\MDn -2.30\ < 0.1GeV/c2. 





neutral B 


charged B [18] 


B{B D^Tr-)B{D^ Dn) 
B{B D^7T-)B{D^ D*7r) 
B{B Dnr-)B{Di D*tt) 


(2.15 ± 0.17 ± 0.29 ± 0.12) x 10"^ 
(2.45 ± 0.42l[J |l[J-39) X 10-4 [19] 
(3.68 ±0.60l°:I™)x 10-4 [19] 


(3.4 ±0.3 ±0.6 ±0.4) X lO'^ 
(1.8 ±0.3 ±0.3 ±0.2) X 10-4 
(6.8 ±0.7 ±1.3 ±0.3) X 10-4 


B{B Dln)B{Dl Dtt) 
B{B D[tt-)B{D[ D*tt) 


(0.60 ± 0.13 ± 0.15 ± 0.22) x 10"^ 
< 0.7 X 10-4 at 90% C.L. [19] 


(6.1 ±0.6 ±0.9 ±1.6) X 10-4 
(5.0 ± 0.4 ± 1.0 ±0.4) X 10-4 



TABLE IV: Comparison of product branching fractions for neutral and charged B decays. 



Taking into account the branching fraction B{f2 — > vrvr) = 0.847lo;oi2 [23] and the corre- 
sponding Clebsch-Gordan coefficients, we obtain 

B{B^ D°f2) = (1.20 ± 0.18 ± 0.21 ± 0.32) x 10-^ 

B{B° D°p°) = (3.19 ± 0.20 ± 0.24 ± 0.38) x 10"^ 

The phases relative to the D2 amphtude are (pp = 2.25 ± 0.19 ± 0.20toll and 4>f^ = 
2.97 ±0.21 ±0.13 ±0.45. 



Results and discussion 



The branching fraction products obtained for the narrow (j = 3/2) resonances are similar 
to the published results for charged B decays as shown in Table IV. The measured values 
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FIG. 8: irn helicity angle distributions for data (points) and MC (histogram). The hatched 
distribution shows the background distribution from the AE sideband region with appropriate 
normahzation. (a) corresponds to the p region jM^r^ — 0.78| < 0.2GeV/c^; (b) the /2 region 
\M^^ - 1.20| < 0.1GeV/c2; (c) the /o region < O.GOGeV/c^. 

of the branching fractions for the broad Dq~^ resonances in neutral B decays are, however, 
significantly lower than those for charged B decays. Preliminary data on 5" D*^tt^'k~ 
decay [19] shown in Table IV indicates a similar behavior for Df and production. One 
possible explanation for this phenomenon is that for charged B decay to D**7r, the amplitude 
receives contributions from both tree and color-suppressed diagrams as shown in Fig. 2. For 
the color-suppressed diagrams, however, D**'s are produced by another mechanism and the 
amplitudes are characterized by the constants /d(3/2) and /d(i/2), with /d(3/2) ^ /d(i/2)- 
The production of the broad resonances and D'^ in charged B decay is amplified by the 
color-suppressed amplitude. As shown in [36], in such a case both ri/2 fits to the sum rule 
and the value of /z)(i/2) are consistent with theoretical estimates. 



CONCLUSION 

A study of neutral 5-meson decays to D^-k^-k^ is reported. We measure the total 
branching fraction of the three-body D^-k^-k^ decays, obtaining B{B^ — >■ D^tt^ti^) = 
(8.4 ± 0.4 ± 0.8) X 10"'^. The intermediate resonant structure of these three-body decays is 
studied. The D^n^n^ final state is described by the production of Dq gvr" with subsequent 
decays ~* Dtt, and also by Dp, -D/2, and Da, where a is a broad scalar (tttt) struc- 
ture. From a Dalitz plot analysis we obtain the mass, width and product of the branching 
fractions for the -D2^: 

M^,+ = (2465.7 ± 1.8 ± 0.8tl^^)MeV /c^, r^*+ = (49.7 ± 3.8 ± 4.1 ± 4.9)MeV, 

B{B^ D;-^tt-) X B{D;+ D\+) = (2.15 ± 0.17 ± 0.29^0.12) x 10^^ 
We observe the production of the broad scalar Dq^ state with the product branching fraction 

B{B° D*Q+7i-) X B{D*Q+ D%+) = (0.60 ± 0.13 ± 0.15 ± 0.22) x 10"^ 

This is the first observation of this decay. The phase of the Dq amplitude relative to that 
of the D2 is determined to be: 

00 = 3.00 ±0.13 ±0.10 ±0.43. 
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The B — > Dp and D/2 branching fractions are measured to be: 

B{B° D^p°) = (3.19 ± 0.20 ± 0.24 ± 0.38) x 10-^ 

B{B^ D°f2) = (1.20 ± 0.18 ± 0.21 ± 0.32) x 10'^, 
and the phases relative to the D2 amphtude are: 

(pp = 2.25 ± 0.19 ± 0.20l[|:^^, 

0/2 = 2.97 ± 0.21 ± 0.13 ± 0.45. 
This is the first observation of the S° — > -D°/2 decay. 
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